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■ A naked singularity occurs in the generic collapse of an inhomogeneous dust ball. We study 

the even-parity mode of gravitational waves from a naked singularity of the Lemaftre-Tolman-Bondi 
spacetime. The wave equations for gravitational waves are solved by numerical integration using the 

' single null coordinate. The result implies that the metric perturbation grows when it approaches the 

Cauchy horizon and diverges there, although the naked singularity is not a strong source of even- 
parity gravitational radiation. Therefore, the Cauchy horizon in this spacetime should be unstable 
^ ' with respect to linear even-parity perturbations. 

m ■ 

: 

I. INTRODUCTION 

The singularity theorems reveal that the occurrence of singularities is a generic property of spacetime in general 
relativity However, these theorems state nothing about the detailed features of the singularities themselves; 

for example, we do not get information from these theorems about whether or not the predicted singularity is naked. 

H ' Here, "naked" means that the singularity is in principle observable. A singularity is a boundary of spacetime. Hence, 
in order to obtain a solution of hyperbolic field equations for matter, gauge fields and spacetime itself in the causal 

bi)' future of a naked singularity, we need to impose a boundary condition on it. However, we do not yet know physically 
reasonable boundary conditions for singularities, and hence to avoid this difficulty, the cosmic censorship hypotheses 
(CCH) proposed by Penrose are often adopted in the analysis of physical phenomena involving strong gravitational 
fields. 

Unfortunately no one has ever succeeded in the proof of any version of the CCH. There is no precise statement of 
CCH which can be readily proved at this time. Given this situation it is worth trying to obtain counterexamples. 
Much effort has been made to search for naked singularity formation in gravitational collapse. 

In the Lemaitre-Tolman-Bondi (LTB) spacetime ^|j7j, a naked shell- focusing singularity appears from generic initial 
data for spherically symmetric configurations of the rest mass density and a specific energy of the dust fluid fs|-|ll|| . The 
initial functions in the most general expandable form have been considered |T^ ] . The matter content in this spacetime 
may satisfy even the dominant energy condition. These results are summarized as follows: In this spacetime, a naked 
singularity appears from generic initial data for spherically symmetric configurations of the rest mass density and a 
specific energy of the dust fluid. Shapiro and Teukolsky numerically studied evolution of collisionless gas spheroids 
with fully general relativistic simulations O]. They found some evidence that prolate spheroids with sufficiently 
elongated initial configurations, and even with some angular momentum, may form naked singularities. Ori and Piran 
numerically examined the structure of self-similar spherical collapse solutions for a perfect fluid with a barotropic 
equation of state 0,^5) . They showed that there is a globally naked singularity in a significant part of the space of 
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self-similar solutions. Joshi and Dwivedi analytically investigated the self-similar spherically symmetric collapse of a 
perfect fluid with a similar equation of sta te |16|] . Harada numerically investigated spherical collapse of a perfect fluid 
without the assumption of self-similarity Jl7|. A spherical cloud of counterrotating particles was investigated by the 
present authors fig] . The spherical gravitational collapse of an imperfect fluid which has only a tangential pressure 
has also been considered ]l9|-p3|]. Further, the naked singularity produced by the gravitational collapse of radiation 
shells |Q and of more general matter were investigated. As for the non-spherically symmetric collapse case, Joshi 
and Krolak revealed that a naked singularity appears also in the Szekeres spacetime with irrotational dust matter 
p6| . The global visibility of this singularity was recently analyzed |2?j . 

In this paper we investigate whether a naked singularity, if such exists, is a strong source of gravitational radiation, 
and we attempt to understand the dynamics and observational meaning of the naked singularity formation. As noted 
above, several researchers have shown that the final fate of gravitational collapse is not always a singularity covered 
by an event horizon. In this case with a small disturbance of spacetime, very short wavelength gravitational waves, 
which are created in the high density region around a singularity, may propagate to the observer outside the dust 
cloud because of the absence of an event horizon. If this is true, extremely high energy phenomena which cannot 
be realized in any high energy experiment on Earth can be observed. Moreover, information regarding the physics 
of so-called 'quantum gravity' may be obtained. Also, these waves may be so intense that they destroy the Cauchy 
horizon. In this paper we consider the generation of gravitational waves during the collapse of a spherical dust ball 
with a small disturbance of the density profile, i.e. perturbations of LTB spacetime. 

Nakamura, Shibata and Nakao |28| have suggested that a naked singularity may emit considerable gravitational wave 
radiation. This was proposed using an estimate of gravitational radiation from a spindle-like naked singularity. They 
modeled the spindle-like naked singularity formation in gravitational collapse using a sequence of general relativistic, 
momentarily static initial data for a prolate spheroid. It should be noted that the system they considered is different 
from that considered in this article and that their result is controversial. There are numerical analyses that may 
support or may not support the results of Nakamura, Shibata and Nakao for prolate collapse [ |l3[ ] and for cylindrical 
collapse HQ. 

Due to the non-linear nature of the problem, it is difficult to analytically solve the Einstein equation. Therefore, 
numerical methods will provide the final tool. However, its singular behavior makes accurate numerical analysis 
very difficult at some stage. In this article, we investigate even-parity linear gravitational waves from the collapse 
of an inhomogeneous spherically symmetric dust cloud. Even for the linearized Einstein equation we must perform 
numerical integration. However, in contrast to the numerical simulation of the full Einstein equation, high precision 
is guaranteed for the numerical integration of the linearized Einstein equation, even in regions with extremely large 
spacetime curvature. Furthermore, the linear stability of known examples of naked singularity formation is necessary 
as a first step to understand the general dynamics near naked singularity formation. 

Recently, Iguchi, Nakao and Harada [^TJ (INH) studied odd-parity metric perturbations around a naked singularity 
in the LTB spacetime. In INH, it was found that the propagation of odd-parity gravitational waves is not affected 
by the collapse of a dust cloud before the formation of the event horizon, even if there appears a central naked 
singularity. The same authors extended their study to consider the generation of gravitational waves from the dust 
collapse including matter perturbation ]3^| . They showed that gauge-invariant variables diverge only at the center, 
and they do not propagate to the outside. For an odd-parity perturbation the evolution of the matter perturbation 
decouples from the evolution of the metric perturbation, while the even-parity matter perturbation couples to the 
metric part. Therefore an even mode seems to be more essential. To investigate the generation of gravitational waves 
in LTB spacetime we should analyze even-parity perturbations. Here we investigate the behavior of the even-parity 
quadrupole metric and matter perturbations in the marginally bound LTB background. We numerically calculate the 
time evolutions of the gauge invariant metric variables. We show that some of metric perturbation variables and the 
Weyl scalar diverge at the Cauchy horizon but that derived the energy flux does not. 

This paper is organized as follows: In Sec. |n| the basic equations are derived; in Sec. Ill the numerical results are 
presented; in Sec. IV we discuss the numerical results; and in Sec. [v] we summarize our results. We adopt geometrized 



units in which c = G = 1 . The signature of the metric tensor and sign convention for the Riemann tensor follow Ref. 
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II. BASIC EQUATIONS 

We consider the evolution of even-parity perturbations of the LTB spacetime to linear order. The background 
LTB spacetime describes the dynamics of an inhomogeneous spherically symmetric dust ball. Using the synchronous 
comoving coordinate system, the line element of the LTB spacetime can be expressed in the form 

ds 2 = g^dx^dx" = -dt 2 + A 2 (t, r)dr 2 + R 2 (t, r)(d9 2 + sin 2 ddcj) 2 ). (2.1) 
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The energy-momentum tensor for the dust fluid is 



T" u = p{t,r)u"u u , (2.2) 

where p(t, r) is the rest mass density and is the 4- velocity of the dust fluid. In the synchronous coordinate system, 
the unit vector field normal to the spacelike hypersurfaces is geodesic, and there is a freedom concerning which timelike 
geodesic field is adopted as the hypersurface unit normal. Using this freedom, we can always set = Sq , since the 
4- velocity of the spherically symmetric dust fluid is tangent to an irrotational timclikc geodesic field. 

Then the Einstein equations and the equation of motion for the dust fluid reduce to the following simple equations: 

A= - = =, (2.3) 

, . 1 1 dF(r) 

R 2 -^ = f(r). (2.5) 
Here f(r) and F(r) are arbitrary functions of the radial c oord inate, r, and the overdot and prime denote partial 



derivatives with respect to t and r, respectively. From Eq. (2.4), F(r) is related to the Misner-Sharp mass function 
|j34j| , m(r), of the dust cloud in the manner 

l-R(t,r) rr pc „\ 

m{r) = 4tt / p(t, r)R 2 dR = 4tt / p(t, r)R 2 R'dr = (2.6) 
Jo Jo 2 

Hence Eq. (|2.5|) might be regarded as the energy equation per unit mass. This means that the other arbitrary function, 
f(r), is recognized as the specific energy of the dust fluid. The motion of the dust cloud is completely specified by the 
function F(r) (or equivalently, the initial distribution of the rest mass density, p) and the specific energy, f(r). When 
we restrict our calculation to the case that the symmetric center, r = 0, is initially regular, the central shell focusing 
singularity is naked if and only if d 2 p\ r= o < is initially satisfied for the marginally bound collapse, f(r) — |p5| , p6[ . 
For collapse that is not marginally bound, there exists a similar condition as an inequality for a value depending on 
the functional forms of F(r) and f(r) 

Next we give a brief introduction to the gauge-invariant formalism of Gerlach and Sengupta (3^j3^] for even- 
parity perturbations around the most general spherically symmetric spacetime. We consider the general spherically 
symmetric spacetime with the metric 

g^dx^dx" = g ab (x d )dx a dx b + R 2 {x d ) lAB {x D )dx A dx B , (2.7) 

and stress-energy tensor 

t^dx»dx u = t ab (x d )dx a dx b + ^t/R 2 (x d )jA B (x D )dx A dx B (2.8) 

where r jABdx A dx B = d6 2 + sin 2 0d<p 2 . Here, lower-case Latin indices refer to radial and time coordinate, while capital 
Latin indices refer to 9 and <fi. 

The even-parity perturbations are 



h ab (x d )Y h a (x d )Y : 



V-( a sym K{x d )R 2 lA a B Y + G(x d )R 2 Z AB i (2 ' <n 



for metric and 



_(At ab (x d )Y At a (x d )Y B \ ( , 

0I ^-[ sym At 3 {x d )R 2 1AB Y + At 2 (x d )Z AB J [Z - W) 

for matter, where Y = Y[ m (x D ) are the scalar spherical harmonics and Z AB — Y :AB + Y-f AB . Here covariant 
derivatives are distinguished as follows: 

lAB-.c = 0, g ab \ c = 0. (2.11) 
For convenience of expression, we introduce 
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Va = R.a/R (2.12) 

and 

Pa = K - \R 2 G^ a . (2.13) 
A set of even-parity gauge-invariant metric perturbations is defined as 

kab = Kb - (Pa\b +Pb\a) (2-14) 

k = K+ " l ' + G - 2v a p a . (2.15) 

A set of even-parity gauge-invariant matter perturbations is defined as 

T a b = At ah - t ab \ c p c - t a c p c \ b - t b c p c \ a , (2.16) 

T a ee At a - t a c Pc - i? 2 (t//4)G, a , (2.17) 

T 3 ee At 3 - (p c /R 2 )(R 2 t//2) iC + 1(1 + l)(t//4)G, (2.18) 

T 2 ee At 2 - (R 2 t//2)G. (2.19) 

The perturbed Einstein equations are expressed only in gauge-invariant perturbations as Eqs. (3.13) of Ref. ]3l| . We 
give these equations in Appendix [A|. 

In this paper we restrict our numerical investigation to the quadrupole mode in the marginally bound background. 
We derive the perturbed equations in th at ca se. Note that, from Eq. (2.3), the background metric variable A is equal 
to R' . Also, we can easily integrate Eq. (2J5) and obtain 

R(t,r)=(—J [io(r)-t] 2/3 , (2.20) 

where to(r) is an arbitrary function of r. The formation time of the naked singularity is to = ^o(O)- Using the freedom 
for the scaling of r, we choose R(Q, r) = r. This scaling of r corresponds to the following choice of to(r): 

t (r) = ^=r 3 ' 2 . (2.21) 

The energy density p is perturbed by adding the scalar term SpY, while the 4-velocity is perturbed by adding 
the term 

8u, t = (V (x d )Y, V 1 {x d )Y, V 2 (x d )Y A ). (2.22) 

The normalization for the 4-velocity yields the relation u^Su^ = 0. This relation implies that Vb vanishes exactly. 
Then there are only three matter perturbation variables, 

T 00 = Sp(t,r), (2.23) 

T 01 =pV 1 (t,r), (2.24) 

T =pV 2 (t,r). (2.25) 

The others exactly vanish: 

Tn =Ti = T 3 = T 2 = 0. (2.26) 

Now we can write down the perturbed Einstein field equations for the background LTB spacetime. The resulting 
linearized Einstein equations are given in Appendix | 



:nd ix |A| . 
5, &~l 



We have obtained seven differential equations, ([A7j)-( A13), for seven variables (four metric and three matter). The 
right-hand sides of four of these equations vanish exactly. Then we can obtain the behavior of the metric variables 
through the integration of them. We transform these equations into more favorable forms. From Eq. (p 



fcoo = -^fcii- (2.27) 
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Using this relation and the remaining equations whose r.h.s. vanish, we obtain evolution equations for gauge-invariant 
metric variables as 



( 2 R"\ , A' . A 



R? H \RR! i?' 3 



R R' 
R~ ~R' 





, , 2R" S R R"R 2R'R' 2R"R' . 
-R" =s — + — =r- + I fcoi 



-R' 



(2.28) 



RR' 



R. R. 
—q - 4—k 
R q R 



RR' 



R' 



R 
R 



(2.29) 



01 



R> k01 



(2.30) 



where q = k — fcoo- If we solve these three equations for some initial data and for the appropriate boundary conditions, 
we c an fol low th e full evolution of the metric perturbations. When we substitute these metric perturbations into Eqs. 
(A7), (A8) and ( A10 ), the matter perturbation variables dp, V% and V2, respectively, are obtained. 

We can also investigate the evolution of the matter perturbations from the linearized conservation equations 
5{T^ V ) = 0. They reduce to 



1 



v 2 



pR 2 R' \ R 



(2.31) 
(2.32) 
(2.33) 



Integration of these equations gives us the time evolution of the matter perturbations. We c an c he ck t he co nsiste ncy 
of the numerical calculation by comparison of these variables and those obtained from Eqs. ( A7), (A8) and (A10). 

To constrain the boundary conditions in our numerical calculation, we should consider the regularity conditions 
at the center. These conditions are obtained from requiring that all tensor quantities be expandable in non-negative 
integer powers of locally Cartesian coordinates near the center ]3^] . The detailed derivation of these conditions is too 
complicated to be presented here. We simply quote the results. The regularity conditions for the metric perturbations 
are 



k~k (t)r 2 , q~q (t)r 4 , k 01 ~k (t)r 3 . 
For the matter perturbations, the regularity conditions at the center are 

Sp~5 Po {t)r 2 , Vx~V 1Q {t)r, V 2 ~V 20 (t)r 2 . 
Therefore all the variables we need to calculate vanish at the center. 



(2.34) 
(2.35) 



III. NUMERICAL METHOD AND RESULTS 



We numerically solved the wave equations ( 2.28 )— ( 2.3C ). Following the method of previous papers, [|3l|,|32| we 
transformed the wave equation (2.28) into the out-going single-null coordinate system. In this section, we present this 
coordinate transformation and explain our background and initial data of the perturbations. In the later half of this 
section, we give our numerical results. 
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A. Numerical method 



In the previous section it was shown that the perturbation variables q, k and fcoi vanish at the center. A careful 
treatment of the differential equations may be required near the center for proper propagation through the center. 
Hence we define the new variables 



qR' 7 /R 4 , k = kR' 4 /R 2 , koi = k 01 R' 5 /R 3 



(3.1) 



These new variables are not identically zero at the regular center and d o not d iverge when they approach the central 
singularity because of the suppression factor R' . We rewri te Eq s. ( 2.28| )-(2.30) in terms of these new variables. 

Next we perform a coordinate transformation for Eq. (2.28) from the synchronous comoving coordinate system 
(t, r) to the single-null coordinate system (it, f), where u is the outgoing null coordinate and f = r. We perform 
the numerical integration of this equation along two characteristic directions. Therefore we use a double null grid 
in the numerical calculation. Whereas we integrate Eqs. (2.29) and (2.30) along the direction r — const. (Detailed 
explanations of the single-null coordinate used in our calculation is given in INH.) As a result, we obtain the first 
order differential equations 



1 d 

— —X = aiX + a 2 W + a 3 Z + a^k + a 5 q, 
a du 



(3.2) 



W = hX + b 2 W + b 3 Z + b 4 k + b 5 q, 



(3.3) 



Z = qX + c 2 W + c 3 Z + ak + c 5 q, 



(3.4) 



k = diX + d 2 W + d 3 Z + d 4 k + d 5 q, 



(3.5) 



df.q = e\X + e 2 W + e 3 Z + e±k + e^q, 
where we have introduced X and W, which are defined by Eqs. (|3.6|) and (3.5), respectively, and 



Z 



R - 



a is given by 



(3.6) 



(3.7) 



a — —. 

u 



The derivatives in Eqs. (3.2) and (3J3) are given by 

d „ df 



du 

dp 



du 



d f = d u 



2R' 



2R' 



-d t + d r = R'dt + d r 



(3.8) 

(3.9) 
(3.10) 



The coefficients a±, a 2 , • ■ • are shown in Appendix |B| Equations (3.2) and (3.6) are integrated along the double-null 
grid. We integrate Eq. (|3.2| ) using the s chem e of an explicit first order difference equation, and we use the trapezoidal 
rule to integrate Eq. ( |3.6|) . Equations (3.3)- (|3.5| ) are integrated along the timelike directions r = const using a first 
order difference method. We interpolate variables to estimate the right-hand sides of Eqs. ( 3.3 )—( 3.5 ) at the same 
radial coordinate r on the previous out-going null slice. 
We adopt the initial rest mass density profile 



P(r) = Pa- 



1 + exp 



1 + exp 



( 2%) 



(3.11) 



G 



where po, r-y and ri are positive constants, and n is a positive even integer. As a result the dust fluid spreads all over 
the space. However, if r 3> fi,r 2 , then p{r) decreases exponentially, so that the dust cloud is divided between the 
dense core region and the envelope, which can be considered as the vacuum region. We define a core radius 

r coie = ri + — . (3-12) 

If we set n = 2, there appears a central naked singularity. This singularity becomes locally or globally naked, 
depending on the parameters po, fi ancl T i- However, if the integer n is greater than 2, the final state of the dust 
cloud is a black hole for all parameter values. Then we consider three different density profiles connected with three 
types of the final state of the dust cloud, globally and locally naked singularities and a black hole. The outgoing null 
coordinate u is chosen so that it agrees with the proper time at the symmetric center. Corresponding parameters are 
given in Table |. Using this density profile, we numerically calculated the total gravitational mass of the dust cloud 
M. In our calculation we adopted the total mass M as the unit of the variables. 
We give the numerical results from the initial conditions for the perturbations 

x= a f q-(e 3 Z + e s q) 
ei 

W = -d 5 q, (3.14) 

it = _gWlWj, (3 . 16) 
64 

«=( 1+ fe) 2 ) 'p* <3 ' 17) 

on the initial null surface. Here k vanishes on this surface and W and Z are diminished near the center. We chose 
r 3 = 0.3r corc . The main results of our numerical investigation do not depend on the detailed choice of the initial 
conditions. 



TABLE I. Parameters of initial density profiles and damped oscillation frequencies, where M = 1. 





final state 


PO 


ri 


ri 


n 


damped oscillation frequency 


(a) 


globally naked 


1 x 10^ 


0.25 


0.5 


2 




(b) 


locally naked 


1 x KT 1 


0.25 


0.5 


2 


0.36+0.096i 


(c) 


black hole 


2 x 1(T 2 


2 


0.4 


4 


0.36+0.093i 
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B. Results 



First we observe the behavior of the metric variables q, k, fcoi and the Weyl scalar, which corresponds to out-going 
waves, 

* 4 = G^ vpa n^fh v n p frf (3.18) 
3 [5 . 2Q k 01 -(k-g)R' 
= -32V^ m ° WW ' (3 - 19) 

where 



2.-2^.0,0) (3.20) 



0,0,-^, = , (3.21) 

outside the dust cloud. The results are plotted in Fig. |l|. We can see that the metric variables q, fcoi and the Weyl 
scalar ^4 diverge when they approach the Cauchy horizon. The asymptotic power indices of these quantities are 
about ~ 0.88. On the other hand the metric quantity k does not diverge when it approaches the Cauchy horizon. The 
energy flux is computed by constructing the Landau-Lifshitz pseudotensor. We can calculate the radiated power of 
gravitational waves from this. The result is given in Appendix For the quadrupole mode, the total radiated power 
becomes 

P=^-k 2 . (3.22) 

The radiated power of the gravitational waves is proportional to the square of k. Therefore the system of spherical 
dust collapse with linear perturbations cannot be expected as a strong source of gravitational waves. 




FIG. 1. Plots of perturbed variables q, k, fcoi and the Weyl scalar ^4 at constant circumferential radius R. The results for 
R — 1, R — 10, R = 100, and R — 200 are plotted. The solid lines represent the results for R — 1, the dotted lines for R = 10, 
the dashed lines for R — 100, and the long dashed lines for R — 200. u = Uo corresponds to the Cauchy horizon. 
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Second we observe the perturbations near the center. The results are plotted in Figs. || and ||. In these figures 
we plot the perturbations at t — io(0) = — 10 _1 , — 10~ 2 , — 1CP 3 , — 1CP 4 , and 0. Before the formation of the naked 
singularity, the perturbations obey the regularity conditions at the center. Each line in these figures displays this 
dependence if the radial coordinate is sufficiently small. In this region, we can also see that all the variables grow 
according to power-laws on the time coordinate along the lines of r = const. The asymptotic behavior of perturbations 
near the central naked singularity is summarized as follows: 

qoc At~ 2 -V, koc At- 1A r 2 , k 01 oc At'^r 3 , 



Sp 
P 



— ocAt^'V, Vi oc At~ 0A r, y 2 ocAr 



where At — to(0) — t. On the time slice at At = 0, perturbations behave as 

goer ' 09 , fcocr" - 74 , k 01 oc r - 92 , 

^cxr" 1 - 4 , Fxocr - 25 , ^ocr 1 - 3 . 
P 

On this slice k and Sp/p diverge and q diverges weakly when they approach the central singularity. On the other 
hand, fcoi and V2 go to zero and V± vanishes slowly. 
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FIG. 2. Plots of perturbed variables q, k and fcoi near the center. The values for At = to — t — 1CP 1 , 10~ 2 , 10 -3 , 10 _4 ,0 are 
plotted. The solid lines represent the results for At = 0, the long dashed lines for At = 10 _1 , the dashed lines for At = 10 -2 , 
the dotted lines for At = 10 -3 , and the dotted dashed lines for At = 10~ 4 . 
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FIG. 3. Plots of perturbed variables Sp,Vi, and Vs near the center. The values for At = 10 _1 , 10 -2 , 10 -3 , 10 -4 , and are 
plotted. The solid lines represent the results for At = 0, the long dashed lines for At — 10 _1 , the dashed lines for At = 10 -2 , 
the dotted lines for At = 10 -3 , and the dotted dashed lines for At = 10 -4 . 



In cases of a locally naked singularity and black hole formation, we expect to observe damped oscillation in the 
asymptotic region outside the dust cloud, as in the odd parity case. The results are plotted in Fig. ||. These figures 
show that damped oscillations are dominant. We read the frequencies and damping rates of these damped oscillations 
from Fig. [| and give them in terms of complex frequencies as 0.36 + 0.096i and 0.36 + 0.093z for locally naked and 
black hole cases, respectively. These results agree well with the fundamental quasi-normal frequency of the quadrupolc 
mode (2Mw = 0.74734+ 0.17792i) @. 

The numerical accuracy of our calculations was checked with the equations that were not used for the derivation 
of Eqs. (2.28|)-(2.30), e.g., Eq. (A10). We define the maximum relative error £ as 



£ = 



ELr)~7 + JK 



( 4 i~ 6 f ) ~k + ^pv 2 



Sjeach term of numerator! 



(3.23) 



We calculated this quantity on the last null surface where the matter variable V% is obtained from the integration of 
Eq. (2.33) using a method similar to that used for Eqs. (3.3)— (3.5). The results are displayed in Fig. @. Exce pt the 
region of small r say, (r < 3 x 10~ 4 ), this value is less than 0.01. Both the numerator and denominator of Eq. ( 3.23| ) 
vanish at the center. Therefore it seems difficult to estimate the numerical errors from Eq. (3.23) when r is small. 
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FIG. 4. Plots of perturbed variables q at constant circumferential radius R — 100 in the locally naked and black hole cases. 
q is normalized with respect to its maximum value, and the origin of the time variable is adjusted to coincide with the time 
when q is maximum. 
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FIG. 5. Maximum relative errors on the last null slice. 



IV. DISCUSSIONS 



In this section we consider the physical interpretation of our numerical results for even-parity perturbations. The 
divergence behavior of the perturbations implies that the linear perturbation analysis near the Cauchy is invalid. This 
fact shows that aspherical effects are important in the naked singularity formation. 
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To consider where these effects are important and what would happen in this region, we should discuss our results 
more carefully. The perturbations grow according to power-laws and diverge only at the Cauchy horizon. Therefore, 
except for the region very near the Cauchy horizon, the perturbations are finite and small when we choose sufficiently 
small initial values. This means that the central region can reach an extremely high density before the breakdown of 
the linear perturbation analysis. While in the region of spacctimc just before the Cauchy horizon, aspherical property 
becomes important for the dynamics of the spacetime. Our results suggest that the Cauchy horizon is unstable and 
that a singularity appears along it. 

The naked singularity of the LTB spacetime is considered as a massless singularity. Gravitational waves, even if 
they have finite energy, would affect the naked singularity. To investigate this effect we should consider back-reaction 
of the gravitational waves. 

For the case of collapse that is not marginally bound, the condition of the appearance of the central naked singularity 
is slightly different from that in the above case (3^]3(| , and hence there is a possibility that the behavior of perturbations 
in this case is different from that in the marginally bound case. However, it is well known that the limiting behavior 
of the metric with t — > to(r) is common to all cases 



Thus we can conjecture that the results of the perturbation analysis for non-marginal collapse would be similar to 
the results for the marginally bound case. 



We have studied the behavior of even-parity perturbations in the LTB spacetime. We have numerically solved the 
linearized Einstein equations for gauge-invariant variables in the case of the quadrupole mode and marginally bound 
background. We have constructed a numerical code which solves the perturbation equations on an out-going single 
null coordinate. For the globally naked case, the perturbed variables q, fcoi and the Weyl scalar grow as powers 
of (uq — u) outside the dust cloud, where the power index is approximately —0.88. Then the Cauchy horizon of this 
spacetime is unstable with respect to linear even-parity perturbations. On the other hand, the perturbed variable k 
is finite just before the crossing of the Cauchy horizon. The energy flux, which is proportional to the square of k, is 
also finite. Therefore inhomogeneous aspherical dust collapse is not expected to be a strong source of gravitational 
wave bursts. 

We have investigated the asymptotic behavior of perturbations near a central naked singularity. If the radial 
coordinate is sufficiently small, the dependence on it is determined by the regularity conditions at the center. Our 
numerical results show this dependence. The time dependence is an inverse power-law in At. At the time of naked 
singularity formation, q, k and Sp/p diverge when they approach the central singularity, while fcoi, Vi an d V? do not. 

For the cases of locally naked and black hole formation, there appear the damped oscillations outside the dust 
cloud. This is consistent with the fundamental quasi-normal frequency of the quadrupole mode of a Schwarzschild 
black hole. 



We would like to thank T. Nakamura for helpful and useful discussions. We are also grateful to H. Sato and 
colleagues in the theoretical astrophysics group at Kyoto University for useful comments and encouragement. This 
work was partly supported by Grants-in-Aid for Scientific Research (No. 9204) and Creative Basic Research (No. 
09NP0801) from the Japanese Ministry of Education, Science, Sport, and Culture. 



The linearized Einstein equations for the spherically symmetric background presented by Gerlach and Sengupta are 




(4.1) 



V. SUMMARY 
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APPENDIX A: LINEARIZED EINSTEIN EQUATIONS 
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-9ab 



1(1 + 1) 1 



R 2 



-9ab 



2k, 



k d + 2 (v a k. b + Vbk, a + fc |b) 

= -167rT Qb , 



k — h I* 



la ^a^c 



i? 2 

-167rT a 



(fc J c + 2t> c fc, c + G/fc) + [kj cld + 2v c kJ d + 2(v c ^ d + v c v d )k cd 



^ab 



k c c l J d + v c k/ lc + nk c c - l V 1 ^k 



= -16ttT 3 , 



k c c = -16ttT 2 , 

where 7Z is the Gaussian curvature of the 2-dimensional submanifold M 2 spanned by x a . Here 

Gab = -2 (v a \ b + v a Vb) + g a b (^vj a + 3v a v a - -^-"j , 

G/ = 2 («J a + v a v a - Tl) . 
For marginally bound LTB spacetime the linearized quadrupole Einstein equations are 



1 



777 <? 



R 

7 q + p9 



R 2 ^ RR' R 
6 (2 R"\ ( R R'\ • 1 , 
'W k+ \RR 1 ~R*) k ~[ 2 R + W) k+ R l2 ~ k 

„ / R RR" R' \ , n R , , 
+2 [WW - RR* + RR 2 fc01 + 2 ^ fc01 = ~ 8 ^' 



^ — Q + { 2^ - 5 I fc' + *' - 4*01 = -87rp7i, 



(Al) 
(A2) 

(A3) 
(A4) 

(A5) 
(A6) 



ir r 



R R' 



R 2 



(A7) 
(A8) 



p/2 p/ p/2 p p/2 p 

„ -ft tt , xl -ft , , Jt it ; .n " 

R 2H R H R H R 

R R jTj j * 



9 + 2 — fc + 2fc + ~ 7572 fc oi' = -lQ7rpV 2 , 



W 3 ' 



1 

R 2 ' 



R' 



R 2 R" _ 
+4 | 



+ fe i + ^7 fc oi = 0, 

R 2 R'\ ; „^2/ „ / RR" it RR'\ , 
+ —) k + 2Rk + 2 \-R^-lV 2 -) km 

_ -R-R , _ / R 2 R" R \ • R 2 , ; 
-2 — -fc i + 2 ^ fc i - 2—kL = 0, 



&oo + 77^ fc n = °- 



Here we have used Eq. ( A13| ) to eliminate ku in Eqs. ( A7)-( Al~2| ). 



(A9) 
(A10) 
(All) 



(A12) 
(A13) 
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APPENDIX B: COEFFICIENTS OF DIFFERENTIAL EQUATIONS 



The coefficients of Eqs. (|3.2[)— (p.q) are 
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RI 3 ^ R 13 



2R< 3 



RR' , RRR' 



R' 4 



i?' 4 ' 



1 

RP' 



= 0, 



rR\C R' 

- 5 R + 6 K- 



1 + R 

rT' 

0, 
1, 
0, 
0, 

r(i + r 



R' 3 

1 

i?' 

0, 

R 1 / 
2 — ^ ( — 



= 0, 



7^-3^(1 
R' R 



R) +5R'. 



(Bl) 
(B2) 
(B3) 
(B4) 

(B5) 
(B6) 

(B7) 
(B8) 

(B9) 
(B10) 

(BH) 
(B12) 

(B13) 
(B14) 

(B15) 
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(B19) 
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APPENDIX C: POWER OF GRAVITATIONAL RADIATION 



In this appendix we calculate the radiated power of the gravitational waves in an attempt to grasp the physical 
meaning of the gauge- invariant quantities M2] . To relate the perturbation of the metric to the radiated gravitational 
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power, it is useful to specialize to the radiation gauge, in which the tetrad components hfg)(6) — W(t>){<t>) an d ^ an 
off as 0(1/R), and all other tetrad components fall off as 0(l/R 2 ) or faster. Note that in vacuum at large distance, 
the spherically symmetric background metric is identical to the Schwarzschild solution, where hereafter we adopt the 
Schwarzschild coordinates, 



ds 2 



1 



2M 
~R 



dr 2 



1 



2M 
~R 



dR 2 + R 2 



+ sin 2 6d(j) 2 ) 



The relation between the line elements Eq. (2.1) and Eq. (CI) is given by the transfer matrix 

1 



dr 



1 - (d t R) 2 
dR = (d t R) r dt + {d r R) t dr. 



{dt+(d r R) t (d t R) r dr}, 



In this gauge, the metric perturbations in Eq. (2.9) behave as 



ha = O ( - 



K = ° { R> 



G = 



g(r-R*) 

R 



° [ R- 2 



where 



R„=R + 2M In 



R 

2M 



1 + const, 



(CI) 

(C2) 
(C3) 

(C4) 
(C5) 
(C6) 
(C7) 

(C8) 



and the out-going wave condition is respected. Then, the gauge-invariant metric perturbations (2.14) and (2.15) are 
calculated as 



k TT = g i2) R + 0(l), 
k Tli = -g^R + 0{l), 
k RR = gWR + 0(l), 



(C9) 
(CIO) 
(Cll) 

(C12) 



where denotes the first derivative of g with respect to its argument. 

In this radiation gauge, the radiated power P per unit solid angle is given by the formula derived by Landau and 
Lifshitz from their stress-energy pseudo-tensor, 



dP _ R 2 
dn ~ 16tt 



dh (0)(4>) \ 1 ( dh {em dh m4>) 
dr J 4 V dr dr 



For the axisymmetric mode, i.e. m = 0, the above formula is reduced to 



where 



MO) 



21 + 1 

47T 



sin 4 9 



d 2 Pi(cost 
(dcos 9)'- 



(C13) 



(C14) 



(C15) 



By using the gauge-invariant quantities and integrating over all solid angles, the formula for the power of the gravi- 
tational radiation is obtained in the following form: 
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(C16) 



(C17) 



where 



(* + 2)! 
(1-2)!- 



(C18) 
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